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I analyze the density perturbations in a cosmological model with a scalar field coupled to ordinary 

matter, such as one obtains in string theory and in conformally transformed scalar-tensor theories. 

- - - ^ The spectrum of multipoles on the last scattering surface and the power spectrum at the present 

Oi ■ are compared with observations to derive bounds on the coupling constant and on the exponential 

y^ ' potential slope. It is found that the acoustic peaks and the power spectrum are strongly sensitive 

^ , , to the model parameters. The models that best fit the galaxy spectrum and satisfy the cluster 

abundance normalization have field energy density Q^ ~ 0.05 — 0.15 and a scale factor expansion 

Ch ■ law a ~ tP, p ~ 0.67 - 0.69. 

^ ; 

1—5 ; 

O ' I. INTRODUCTION 

^ ' Perhaps the most important concept in modern cosmology is that fundamental physics, along with gravity, shapes 

C^ , the distribution of matter on very large scales. Fundamental physics enters in at least two distinct ways: through the 
t^ ■ potential of the inflationary field, which sets the initial conditions of the fluctuation field, and through the properties 
of the dark matter, which govern the evolution of the fluctuations up to the present. As a consequence, the imprint 
of the density fluctuations on the microwave background and on the galaxy distribution allows tests of basic laws of 
^^ ' physics that, in many cases, could not be realized with any other mean. 
^N^ . An impressive array of different proposals have been formulated for as concerns the inflationary side of the story, that 

is, the initial conditions. So far, there is not an overwhelming reason to modify the simplest inflationary prescription 
of a flat spectrum, although several variations on the theme, like a small tilting (Lucchin & Matarrese 1985, Cen et 



o 
o 



I ' al. 1992) or some break in the scale invariance (Gottloeber et al. 1994, Amendola et al. 1995) or the contribute of 
O I primordial voids (Amendola et al. 1996) cannot be excluded either. 

-4-^ ■ Similarly, many theories have been proposed for as concerns the evolution of the fluctuations, trying to elucidate 

^ ] the nature and properties of the dark matter component. A partial list of the dark matter recipes includes the 

• • . standard CDM and variations such as CDM plus a hot component (MDM) , or plus a cosmological constant ( ACDM) , 

_ ^ or plus a scalar field (0 CDM). The latter class of models, in particular, has been explored greatly in recent times, to 

k^ „ various purposes. First, a light scalar field is predicted by many fundamental theories (string theory, pseudo-Nambu- 

^ ■ Goldstone model, Brans-Dicke theory etc), so that it is natural to look at its cosmological consequences (Wetterich 

5r ' 1995, Frieman et al. 1995, Ferreira & Joyce 1998). Second, a scalar field may produce an effective cosmological 

constant, with the benefit that its dynamics can be linked to some underlying theory, or can help escape the strong 

constraints on a true cosmological constant (Coble, Dodelson & Frieman 1997, Waga & Miceli 1999). In turn, this 

effective cosmological constant can be tuned to explain the observation of an accelerated expansion (Perlmutter et 

al. 1998, Riess et al. 1998) and to fix the standard CDM spectrum as well (Zlatev et al. 1998, Caldwell et al. 1998, 

Perrotta & Baccigalupi 1999, Viana & Liddle 1998). Finally, even a small amount of scalar field density may give a 

detectable contribution to the standard CDM scenario, similar to what one has in the MDM model (Ferreira & Joyce 

1998, hereinafter FJ). 

In this paper we pursue the investigation of the effects of a scalar field in cosmology by adding a coupling between 
the field and ordinary matter. Such a coupling has been proposed and studied several times in the past (e.g. Ellis et 
al. 1989, Wetterich 1995) but, as far as we know, its consequences on the microwave background and on the power 
spectrum have not been determined. Scope of this paper is to solve the fluctuation equations for a coupled scalar 
field theory, and to compare with the already available data in the microwave sky and in galaxy surveys. We refer to 
this model as coupled (j)CDM. Up to a conformal transformation, the model we study is equivalent to a non-minimal 
coupling theory in which the scalar field couples to gravity like in a Brans-Dicke Lagrangian; the perturbations in 
such models have been investigated by Chen & Kamionkowsky (1999) and Baccigalupi et al. (1999) in a background 
in which the scalar field acts like a dynamical cosmological constant (see also Uzan 1999). The model we present here 
is in fact more general, since a wide class of non-minimal coupling models can be recast in the form we study below 
(Amendola et al. 1993, Wetterich 1995, Amendola 1999). 



There are several models of cosmological scalar field in the literature, essentially characterized by the scalar field 
potential and by the initial conditions. We can divide the models in two broad class: in the first one, the field potential 
energy dominates at the present, so that it resembles a cosmological constant. In the second one, the field kinetic 
energy is not negligible, and the field adds to the ordinary matter as an additional component, like in MDM models. 
To this second class belongs the model of FJ. They adopt an exponential potential for the scalar field, able to drive an 
attractor scaling solution which self-adjusts to the dominant matter component. In such a model, the density fraction 
of the field does not depends on the initial conditions, but is determined by the potential parameters. Therefore, the 
coincidence that the energy density in the field and in the matter components are comparable can be explained by the 
underlying physics (the field potential) rather than by the initial conditions. Although the coupling we will introduce 
can be applied to any scalar field model, we focus our attention in this paper on the exponential potential model of 
FJ. Beside being particularly simple, because of its attractor properties (Wetterich 1988, Ratra & Peebles 1988), such 
a model is also easily falsifiable, because the effect of the scalar field is important at all times (not just at the present 
as when the field acts as a cosmological constant), and therefore induces a strong effect on the cosmological sky. As 
a consequence, the constraints we derive on the model parameters are rather strong. 

The same exponential potential also allows solutions which belong to the first class mentioned above, in which the 
field acts much like a cosmological constant, and drives at the present an accelerated expansion. These solutions, 
and their linear perturbations, have been studied by Viana & Liddle (1998) and Caldwell et al. (1998). The effect of 
adding a coupling to these models will be analyzed in a subsequent work. 

II. COUPLED SCALAR FIELD MODEL 

Consider two components, a scalar field (j) and ordinary matter (e.g., baryons plus CDM) described by the energy- 
momentum tensors T^^i^\ and T^^i„-^y General covariance requires the conservation of their sum, so that it is possible 
to consider a coupling such that, for instance, 

T^^m);, = -CT^m)<^;,. (1) 

Such a coupling arises for instance in string theory, or after a conformal transformation of Brans-Dicke theory (Wet- 
terich 1995, Amendola 1999). It has also been proposed to explain 'fifth-force' experiments, since it corresponds to 
a new interaction that can compete with gravity and be material-dependent. The coupling arises from Lagrangian 
terms of the form (Wetterich 1995) 

- ml exp(-CK(/))V',pV''^, (2) 

where k^ — 8ttG and ip is the ordinary matter field of mass m^, e.g., the nucleon field. 

The specific coupling (m) is only one of the possible form. Non-linear couplings as CTi^^\F{(f))(f);i^i_ or more complicate 
functions are also possible. Also, one can think of different coupling to different matter species, for instance coupling 
the scalar field only to dark matter and not to baryons. Such a species-dependent coupling has been proposed by 
Damour, Gibbons & Gundlach (1990), and shown to be observationally viable. Notice that the couphng to radiation 
(subscript 7) vanishes, since T(.^) = 0. Here we restrict ourselves to the simplest possibility (^, which is also the 
same investigated earlier by Wetterich (1995) and is the kind of coupling that arises from Brans-Dicke models. For 
instance, a field with coupling to gravity ^S^cjpR in the Lagrangian acquires, after conformal transformation, a coupling 
to matter of the form (nl) with 

^= ^^77; (3) 

(6 + 1/0'/' 

in the limit of small positive coupling this becomes 

C = ^yf. (4) 

There are several constraints on the coupling constant C along with constraints on the mass of the scalar field 
particles, reviewed by Ellis et al. (1989). The constraints arise from a variety of observations, ranging from Cavendish- 
type experiments, to primordial nucleosynthesis bounds, to stellar structure, etc etc. Most of them, however, apply 
only if the scalar couples to baryons, which is not necessarily the case, and/or involve the mass of the scalar field 
particles, which is unknown. The most stringent bound, quoted by Wetterich (1995) amounts to 



\C\ < O.lMpi (5) 

but again holds only for a coupling to baryons. Moreover, these constraints are local both in space and time, and could 
be easily escaped by a time-dependent coupling constant. In the following we leave therefore C as a free parameter. 
The constraints from nucleosynthesis refer to the energy density in the scalar component. This has to be small 
enough not to perturb element production, so that at the epoch of nucleosynthesis (Wettcrich 1995, Sarkar 1996, FJ) 

f^0<O.l-O.2 (6) 

We will see that this bound is satisfied by all the interesting models. 

There is an immediate consequence of the coupling for as concerns cosmology. The coupling modifies the conser- 
vation equation for the ordinary matter, leading to a different effective equation of state for the matter. This alters 
the scale factor expansion law in matter dominated era (MDE) from a ~ t^^^ to t^, p ^ 2/3 . In turn, this has 
three effects. First, the sound horizon at decoupling (when decoupling occurs in MDE) is modified with respect to 
the uncoupled case, being larger for p > 2/3 and smaller in the opposite case, as we will show. This modifies the 
peak structure of the microwave background multipolcs. Second, the epoch of matter/radiation equivalence moves to 
a later epoch if p > 2/3 and to an earlier epoch in the opposite case. This shifts the range of scales for which there 
is the growth suppression of the sub-horizon modes in radiation dominated era (RDE), leading to a turnaround of 
the power spectrum on larger scales ii p > 2/3 (smaller if p < 2/3). Finally, the different scale factor law modifies 
the fluctuation growth for sub-horizon modes in the MDE, generally reducing the growth for all values of C. Similar 
effects have been observed by Chen & Kamionkowsky (1999) and Baccigalupi et al. (1999) in Brans-Dicke models. 
The next sections investigate these effects in detail. 

III. BACKGROUND 

Here we derive the background equations in the conformal FRW metric 

ds^ = a^{-dT^ + S^jdx'dx^). (7) 

The scalar field equation is 

4> + 2H<i) + a^U,^^Cpma^, (8) 

where H = d/a , and we adopt the exponential potential 

Uicj)) = Ae"l'. (9) 

The matter (subscript to) and the radiation (subscript 7) equations are 

Pm + "iHpm = -Cpm4' (10) 

Py + AHp^^O. (11) 

Denoting with tq the conformal time today, let us put 

a(To) = 1, Pm(To) = -^f^m = PmO, Pjiro) =^ p^o, 0(to) = 0o- (12) 

an 

Without loss of generality, the scalar field can be rescaled by a constant quantity, by a suitable redefinition of the 
potential constant A. We put then ^o = 0. This gives 

Pm = Pmoa~^e~"^'^, 
Pi = Pioa^"^- (13) 

The (0,0) Einstein equation can be written 



2 _ K f PmO -c<t> I P-rO 1 



i,2 , rr 2 



3 \ a o'^ 2 J 

The dynamics of the model is very simple to study in the regime in which cither matter or radiation dominates. 
Assume that the dominant component has equation of state 



p^ {w~ l)p. 



(15) 



Then, following Copeland et al. (1997) we define 



V6H' 



y = 



kVu 



(16) 



and introduce the independent variable a = log a. Notice that x'^ and y^ give the fraction of total energy density 
carried by the scalar field kinetic and potential energy, respectively . Then, we can rewrite the equations as (Amendola 
1999) 



x' = —3x + 3x 
y' = fixy + 3y 



2,1/1 2 2\ 

X + -w{l -X -y ) 

2 I 1 /-I 2 2\ 

X + -w[l ~x ~y ) 



-fiy^ + f3il~x^-y^), 



where the prime denotes here d/da and where we introduce the adimensional constants 



P- 



3C 

2"^' 



M 



3s 

2k' 



(17) 



(18) 



(in Amendola 1999 we defined /3 one half of the definition above). Notice that in this simplified system with a single 
component, plus the scalar field, the constant (3 is the coupling constant for the dominant component only, so that 
we are implicitly assuming (3 = during RDE. The system is invariant under the change of sign of y and of a. Since 
it is also limited by the condition p > to the circle a;^ + y^ < 1, we limit the analysis only to the unitary semicircle 



of positive y. The critical points, 
field equation of state is 



i.e. the points that verify x' 



y 



0, are scaling solutions, on which the scalar 



2x' 



w^ 



x^ + y^ 



const, 



(19) 



the scalar field total energy density is fi^ = x^ + y^ , and the scale factor is 



= tP, 



^^3 



^4,{w^ ■ 



(20) 



(t being the time defined by dt = a{T)dT). 

The system ( fij\ ) with an exponential potential has up to five critical points, that can be classified according to the 
dominant energy density: one dominated by the scalar field total energy density, one in which the fractions of energy 
density in the matter and in the field are both non-zero, one in which the matter field and the field kinetic energy 
are both non-zero, while the field potential energy vanishes, and finally two dominated by the kinetic energy of the 
scalar field (at 2; = ±1) . The critical points are hsted in Tab. I, where we put g{(3, w, fi) = 4/3^ -I- 4/3/x -t- 18w. For any 
value of the parameters there is one and only one stable critical point (attractor). More details on the phase space 
dynamics in Amendola (1999) and, for /9 = 0, in Copeland et al. (1997). 
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Tab. I 



The perturbations on solutions converging toward the attractor a have been studied in Viana & Liddle (1998) and 
in Caldwell et al. (1999) for zero coupling. In this case the scalar field is starting to dominate today, and mimics a 
cosmological constant. The case of interest here is instead the solution h in Tab. I, since this is the only critical point 
which allow a partition of the energy between the scalar field and the matter and (contrary to c) is stable also in the 



RDE, when /? = 0. The solution b is compatible with a p larger or smaller than 2/3. It exists and is stable (that is, 
is an attractor) in the region delimited by fi < /i_ and fi > /i+ and the two branches of the curve 

Mo = -^(4/32 + 18«;-W). (21) 

The scale factor slope on the attractor is (Wetterich 1995) 

and, ii w = 1, is inflationary for 

2/3 > fi. (23) 

The parametric space region in which the attractor exists is shown in Fig. 1. For any value of the parameters /3,/i 
there is a pair of observables 51^, p. When radiation dominates, P — 0, and the scale factor is the usual RDE one, 
p = 1/2. The mapping from /3,/i to fl^,p is shown in the same Fig. 1: as one can see, to get a large p a large fJ^ 
is also needed. In Fig. 2 we show the phase space of the system for Q^ = 0.1, p — 0.7 assuming matter domination. 
Notice that only for /3 y^ there is the possibility to get an inflationary attractor with Q^ < 1, as some observations 
suggest. It can be easily demonstrated that the coupled exponential potential with 2/3 > /x is the only model that 
allows inflationary attractors with a non-vanishing matter component. Although such a possibility is intriguing, it 
is hardly realistic, since an inflationary expansion that lasted for most of the MDE would not allow any fluctuation 
growth via gravitational instability. 

When both radiation and matter are present, the system goes rapidly from the radiation attractor, for which 

n^iR)^6/^J.^ PR = 1/2, (24) 

to the matter attractor 

9 _2 f P 



^4>= ... ..,2 ^ P-o(l + -)- (25) 



It is convenient to note that /3//i is a measure of the deviation from the uncoupled po = 2/3 law in MDE: 

^ = ^ = ^-1 = ^-1 = ^. (26) 

/i s 2 Pa p 

We give also the relation between the parameters (/3,/i) and the observables (rj0,p) : 

18pri^ 

T \ 1/2 1 r / „ \ 1/2 1 

(27) 
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small ^ we have 
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C-KVSf^T^/'^. (28) 
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Since the slope and the matter content in the MDE depend on the model parameters, the equivalence epoch 
(subscript e) also depends on them. It is easy to see that the following relation holds 

^4-3po/p = fli, (29) 

PmO 

Clearly, the equivalence occurs earlier with respect to the uncoupled case if p < 2/3 (that is, C/s < 0), later Up > 2/3 
(that is, C/s > 0). 

We will make often use of the fact that on the attractor in the RDE (subscript R) and in the MDE (subscript M) 
we have 



= aM,R log a, 



(30) 



where 



Finally, it is useful to note that 



4 3 

s s + C 



CoiM — : — :::; — ~3 — 



s + C p 



(31) 



(32) 



(the latter is valid for ^ < 1). 



IV. PERTURBATIONS 

We now proceed to study the evolution of the perturbations in the coupled 0CDM theory. This involves the 
following tasks: 1) calculate the linear perturbation equations (we choose the synchronous gauge for the perturbed 
metric) for the coupled system of baryons (subscript fe), CDM (c), radiation (7), scalar field (0), massless neutrinos 
(v); 2) establish initial conditions (we adopt adiabatic initial conditions); 3) evolve the equations from deep into the 
radiation era and outside the horizon down to the present; 4) calculate the radiation fluctuations on the microwave 
background and the matter power spectrum at the present; 5) compare with observations. 

Let us identify the effects of adding a scalar field to standard CDM. The field component clearly induces two main 
consequences for as concerns the perturbation equations: delays the epoch of equivalence, because the matter density 
at the present is smaller than without scalar field, and changes the perturbation equations. The first effect induces 
a turn-over of the power spectrum at larger scales, just as in the case of an open universe, or a model with a large 
cosmological constant, so that the power spectrum normalized to COBE has less power on small scales, as observed. 
The modification to the perturbation equations goes in the same direction: the evolution in the MDE for sub-horizon 
modes is suppressed with respect to standard CDM, as we will see below. The evolution equations in the other 
cases (super-horizon modes, RDE) give the same behavior as for the pure CDM . The net result is that FJ find that 
fi^ ~ 0.1 gives a good fit to observations, comparable or superior to MDM or ACDM. 

When we insert the coupling, the two effects above mentioned are again the dominant ones. But now, the conse- 
quences of the coupling can be in either directions, that is, the equivalence epoch can be delayed or anticipated, and 
the perturbations can be either suppressed or enhanced with respect to the uncoupled case, although not by a large 
factor. To understand this effects we first discuss analytically the perturbation equations. Following the discussion 
in FJ, we simplify the problem by reducing the system to three components, CDM, scalar field, and radiation. The 
notation is 



6 = 5p/ p, ip — 5(j), Vi — Sui, ik^Vi 



(33) 



where Ui is the comoving velocity. The perturbation equations in synchronous gauge are: 
Scalar field equation: 

<P + 2Hip + k^p + a^U_^(f,p + —h4> = Ca^pm^cSc, 



(34) 



CDM: 



6c = -de ^ ^^ ^ C''/'' 
e^ = -HOc + C{k^p - 



(35) 
(36) 



Radiation: 



4 2- 

5^ = --Oj - -h, 

. kK 



(37) 
(38) 



Energy- momentum tensor: 



a'^ST^ = a-" {ScPc + S^Pj) + c/}^ + a^U,ip (39) 

2 2 

^tk'ST^ = ^w-,d^p^ + 0^ (40) 

,2rT^J _ „2r „ o I 1,- „2 



a^(57^» = -a^J^p^ - 3 U(^ - a t^,^ ) ■ (41) 
Metric: 

nil = 2k^T] + Sna^STl^ (42) 

7) = 47r^7fc'<5rO (43) 

h = -Hh-8Tra^{ST^ ~6T,'). (44) 
Deriving Eq. (^5|) and inserting equation (Q) we get 

3 

—J 
2 



^c + -ff(5c - -H^injc + 2n^S-,) - 8tt ( 2<j>Lp - sa^Uip) + C(Hlp + k^Lp + Cp- AH^n^if) = 0. (45) 



The equation for the scalar field becomes (putting 9c = 0) 



/o rr2 

tf + 2Hip + k^if + s^a^Uip ^Sc^^C { ^—ricSc + ^'p] ■ (46) 



Finally, the radiation equation is 

6^ + y ^7 - ^ (^c + 0c + C^) = 0. (47) 

The adiabatic initial condition gives now, putting for the initial value of the scalar field ip = x^c (x is determined 
below) , 

S^^P^il + Cx). (48) 

In the large scale limit, fc^ — > 0, and in RDE, where H — r^^ and flc -^ 0, and assuming the adiabatic condition, the 
system reduces to 

^c + T^^Sc - 4T-^Scfl-y - Svr h^ifi - sa^Utp) + C{T~^ip + (p ~ At-^Q.^lp) ^ 0, 

Cp + 2T-^Lp + s^a'^U^-5c^-C^Lp = Q. (49) 

Inserting the RDE attractor solution for 0, we obtain that the growing mode both for 5 and ip goes as r^ . Therefore, 
the super-horizon perturbations in RDE grow similarly in CDM, in (/)CDM, and in coupled 0CDM. Moreover, we have 
that, initially, 

4 / 4C^ "^ 



^ = --^. 1 + -) ^x5.. (50) 



Therefore, the initial condition for the CDM density fluctuations on the attractor in the RDE is 

1 h 



21 + Cx 
Now we consider the super- horizon modes in MDE. The equations are now 

^ '72o X o_ (r,),- „„2, 



(51) 



5c + H5c - -H-'flcSc - Stt ( 2(t)tp - sa^Up ) + C{Hlp + ^) = 0, 

'31/2 



^ + 2Hip + s^a^Up - (5c0 - C I ——VtcSc + 4)^] =0. (52) 

The growing mode is again r^, that is, there is no difference with respect to the standard case. 



In the sub-horizon regime, neglecting the gravitational feed-back, we have in RDE 



Sc + H6c - 4H^n^S^ - Stt (2,j)ip - sa^Uip) + C{Hip + k^Lp + <p~ AH^Q^cp) = 



^ + H{2-Ca)ip + k^ip = 



6j + —6j — 0. 
o 



The oscillating behavior of ip and of S-y gives a negligible influence on 6c-, so that 

L + hL = 0. 



(53) 
(54) 

(55) 
(56) 



which gives 5c = const, logr, once again with no difference with respect to standard CDM. 

We finally come to the regime where the new physics makes the difference. In the sub-horizon MDE regime, 
neglecting again the gravitational feed-back, we have 



-H'^VlcSc - Svr 2(/)(^ - sa^Uip 



H5c 



Neglecting the oscillating behavior of p we obtain 



C 



-Hip + C [ -^^—flcSc + (fxp 







(p + 2Hip + k^p- C 



3H' 



ficSc + 4"p I = 0. 



Sc + H{1 + CaM)5c ~ ^H^nc5c{l - ^) = 0. 

Inserting the trial solution &c — Ba"^ we obtain two solutions for m: 

1-p 



m± 
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[-i±[i + F{n^,p)f'} 



where 



F{n4„p) 



6p{l - n^) [-8 + 26p + 3in^ - 7)p2] 



(p-l)2[2 + 3p(17^-l)] 



For p — 2/3 this reduces to the form found in FJ 



m±-i(-l±y25" 



2m, 



(57) 
(58) 

(59) 
(60) 
(61) 
(62) 



In Fig. 3 we show the contour plot of m^{fl^,p). This figure is crucial for the understanding of the perturbation 
evolution, so we discuss it at some length. First, we observe that for all values of i~i^,p there is suppression with 
respect to CDM: the slope is in fact always less than 1 and, for p ~ 2/3, the slope is smaller for larger fi^. Second, 
we notice the unexpected fact that the value p = 2/3 is close to the maximum for all values of 17^, and closest for 
small n^. For fi^ = 0.1, for instance, the maximum is at p = 0.672, while for fi^ — 0.6 it is at p = 0.710. This 
implies immediately that the coupling does not enhance much the fluctuation growth with respect to the uncoupled 
case, while it can sensibly reduce it further as long as p is far from 2/3. Third, there is only a finite range of p, almost 
centered around 2/3, for which real values of m± exist. Beyond that range, the power-law solutions of Eq. (|5^) are 
replaced by oscillating solutions cos (logo), in which the restoring force is the coupling interaction. 

Let us then summarize the asymptotic evolution of the fluctuations in the coupled model. There are two relevant 
cases. If p > 2/3, the equivalence epoch occurs later than in the uncoupled case. Then, smaller wavenumbers reenter 
during the RDE than in the uncoupled case, and therefore there is extra suppression at these scales. Then, in the 
subsequent MDE regime, the modes are further suppressed with respect to the uncoupled case, unless p is close to 
2/3. The transfer function will be then more steeply declining with respect to the uncoupled case. If p < 2/3, on the 
other hand, the equivalence occurs earlier, and the scales smaller than 27rre/ae are less suppressed. At the same time, 
the MDE regime induces again a slower fluctuation growth, so that there is an intermediate region of wavenumbers 
with a depleted transfer function, and a large wavenumber region with an enhanced transfer function. Fig. 4 displays 
some of these features. 



The only important difference that arises when the baryons are added is in the tight couphng approximation. 
Referring to the notation used in Ma & Bertschinger (1995), we have the two equations 



, . , 1 . \ 1 



e^ = k^ i-s^ - C7^ \ + — [Ob - e^) (63) 

Bt, = -HOi, + clPSb --{9b- e.,) + C(A:V + ^^fe)- (64) 

Tc 

The shp equation 6^ — 9^ in the tight couphng approximation can be derived exactly as in Ma & Bertschinger (1995), 
taking into account that now (here n^ is the electron density and (Jt the Thomson cross section) 



{an^cjTY^ - fl-^eC^-^, f^ = {2H + cA Tc 



R=^, R=(h-C<J))r. (65) 

To second order in Tc we obtain that the slip between baryons and photons is 
9b -9^ = {9b - 9y) 

c(HPip + H^9b + P(p + 4>9b)]. (66) 

The equation for the photons is 

'1 



-R-' 



9b + H9b-k'ci5b-Cik^ip + (l)9b) + k' [ -S^ - a^ ] . (67) 



This concludes the analysis of the asymptotic regimes in the coupled </> CDM model. The results that will be 
presented in the next Sections make use of the full machinery of the Boltzmann code, as implemented in the CMBFAST 
code of Seljak and Zaldarriaga (1996), opportunely modified to take into account the coupled scalar field (including 
the transient from the RDE attractor to the MDE one). The equations are essentially the same as in FJ, with the 
new terms due to the coupling as detailed above. We tested the code with the results of FJ when p = 2/3, and we 
also checked our results with the asymptotics found above. 

V. COMPARISON WITH OBSERVATIONS: COSMIC MICROWAVE BACKGROUND 

The main effect of the coupling on the cosmic microwave background is on the location and amplitude of the 
acoustic peaks. The location of the peak is related to the size of the sound horizon at decoupling (subscript d). Since 
the photon-baryon fluid has sound velocity 



1 R 

— T r — 

3 ' ~ w'+R' 



c%b^^r, r ^ f\ ^ , (68) 



where w'^ ^ I + Ca/'i, the sound horizon is 

dr r" da 



(3r)^/^ Ja {iH^a-^rf^' 



(69) 



This expression can be simplified as follows. First, we put ourselves in the case Ce ^ a^ — 10 ^ and neglect the RDE 
stage altogether. In MDE we have 

H^a^ ~ Hlae-^'t'. (70) 

Then we can write, remembering that on the attractor e^*^*^ = a^*^", and defining the standard sound horizon 
r,o = 2ay'i7o~V\/3 



/o 2(ai-C'"r) 
We can further simplify, for r ~ 1, i.e. i? 3> 1 (which is true at decouphng) 



rso /■"'' da , , 

" ' (71) 



Cq/2 



and the corresponding peak multipole is, for (5p/p <C 1 



«,,,fe c. -^ ^ 4 (1 + Ca) a;^"/' = 4 (l - 3^) a^'^^''/^ (73) 



where the standard peak multipole is 



4 = —rr - 200. (74) 



The qualitative behavior is clear: for p < 2/3 there is a larger t^^ak than in the uncoupled model, for p > 2/3 a 
smaller (.peak- For instance, for p = 0.65 we expect (.peak — 250, in agreement with the numerical results. 

We calculated the Ci spectrum for several coupled (/)CDM models, parametrized by the two observables f^^jp. The 
range of values we explore, in this and in the next Section, is 

17^ = 0.05-0.2, p = 0.65 -0.70. (75) 

The values of the other relevant parameters are fixed as follows 

h = 0.7, VLi, = 0.04, f^A = 0, n = 1. (76) 

In Fig. 5 we display the multipole spectra. As anticipated, the acoustic peaks move to larger multipoles as p decreases. 

There are two other effects worth discussing: the amplitude of the acoustic peaks and the slope of the multipole 
spectrum at small (. The amplitude of the peak is depressed as 5p/p increases, save for values close to 2/3, since the 
matter fluctuations that drive the radiation peaks are suppressed, as shown above. The small ( region is dominated 
by the Sachs- Wolfe (SW) effect. As well known, the integrated SW (ISW) effect in flat space vanishes only if the 
fluctuations grow as a, which is not the case here. The ISW then adds at small multipoles and tilts the Ci spectrum. 
Moreover, the overall normalization now takes into account the ISW power, and as a consequence the normalization 
for the perturbation at decoupling time is reduced. This effect shows also in the final amplitude of the power spectrum. 

Deriving precise constraints from the whole set of observations on the CMB requires considerable detail in the 
statistical procedure, beyond the scope of this paper. Here we content ourselves to derive rough limits on the 
parameters. It is probably safe to state that current observations rule out values p < 0.63 oi p > 0.72, although the 
present level of errors does not permit to attach a strong significance to such bounds. Future precision observations 
around the first peaks are likely to constrain p to two decimal digits. As already found by FJ, on the other hand, 
the microwave sky does not impose strong constraints on Q,^, since this parameter infiuences mainly the fiuctuation 
growth, and thus the absolute normalization. To constrain it, we have to evaluate the present power spectrum of the 
fluctuations. 

VI. COMPARISON WITH OBSERVATIONS: POWER SPECTRUM 

The analytical expression ( |6l| ) for the fluctuation growth exponent found in Section 3 is a clear guide to the results 
of this Section. As anticipated, the coupling introduces an extra suppression for the scales that enter the horizon in 
the MDE. The suppression is larger for models with high Q^ and high | Sp/p \ . A small suppression factor, as well 
known, helps to bring the standard CDM model into agreement with observations. FJ found that the best uncoupled 
models have il^ — 0.1; here we see that the coupling allows also models with smaller il^, but Sp/p ^ 0, to meet the 
observations. This can be helpful to reduce the constraints from nucleosynthesis, which, in some restrictive analysis, 
require fi^ < 0.1. 

In Fig. 6 we report the power spectra A'^{k) — k'^ P (k) / (/Itt^) normalized to COBE, compared to the data as 
compiled and corrected for redshift and non-linear distortions by Peacock & Dodds (1994). For a quantitative 
comparison, we plot in Fig. 7 the contour plot of a^{Q,fj,,p), the number density variance in 8 Mpc/h. spheres. The 
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models with a value of erg larger than 0.5, as required by cluster abundance (White, Efstathiou & Frenk 1993, Viana 
& Liddle 1996, Girardi et al. 1998), have fl^ < 0.15 and at most a small deviations from p ^ 2/3. The suppression of 
as with respect to COBE-normalized standard CDM is due both to the growth suppression in MDE and to the fact 
that now the COBE normalization includes the ISW effect. 

For as concerns the shape of the spectrum, the comparison with the galaxy data is uncertain due to biasing. 
Assuming a scale-independent bias between matter and galaxies, we can quantify the agreement with the data by 
evaluating the x^ of the ratio between the theoretical and the galaxy spectrum, Ri = PG{ki)/ Pxiki), that is by 
evaluating 

X^ = Y^[r,-R^ /a\h), (77) 

i 

where R, a'^{k) are the average and variance of Ri, neglecting cosmic variance. The contour plots of x^{^4>iP) ^^^ in 
Fig. 8. They show, as anticipated, that the models with p > 2/3 follow better the real data because are more bent 
at small scales. The best models among those studied here have x^ — 13 for N = 15 degrees of freedom (16 real 
data, minus the average R estimated from the data themselves). For instance the model with il^ — 0.15,p = 0.685 
gives a very good fit, and has erg — 0.6 as required. Notice that we performed the fits without varying all the other 
cosmological parameters, which, at least in principle, can be determined by other observations. 

In Fig. 9 we summarize the constraints obtained in this Section, considering the models which have erg G (0.45, 0.75) 
and x^/A^ < 2. The cluster abundance normahzation usually quoted for O^ ~ 1 is erg ~ 0.60 ± 0.05, but this 
is calculated for standard models, so that conservatively a larger region has been adopted. Only a stripe around 
p — 0.67 — 0.69 and fi^ = 0.05 — 0.2 appears to be viable, if the bias is indeed scale-independent. From Eq. ( [2^ ) we 
deduce a limit 

< C < IM-i (78) 

which, although still far from the limit O.lMp^ quoted in Wetterich (1995) from local measures, is global and applies 
even if the scalar field is not coupled to baryons, as proposed in Damour, Gibbons & Gundlach (1990). Future data 
will certainly tighten the constraint even further. 

VII. CONCLUSION 

In this paper we discussed the perturbations of a coupled scalar field with exponential potential on the CMB and 
on the present large scale structure along an attractor solution. For as concerns the CMB, we found that the coupling 
induces a strong effect on the location and amplitude of the acoustic peaks, due to the variation of the scale factor 
expansion law. Future precision measures in the region £ > 100 have the potential of constraining the coupling to two 
orders of magnitudes better than at the present (see for instance the discussion in Chen & Kamionkowsky 1999). 

We found that subhorizon perturbations are always more suppressed in MDE with respect to standard CDM, no 
matter what the parameters fi^ and p are. Moreover, the suppression increases for p far from the standard value. The 
amplitude erg at the present is between 0.5 and 1 only for fJ^ e (0.05 — 0.2), being smaller for larger fi^, as already 
found in the uncoupled case by FJ. Adding the request to fit the galaxy power spectrum shape, the parametric space 
is reduced as in Fig. 9. A positive coupling has the advantage to warp the spectrum to a closer agreement with the 
data. 

The background solution we adopted here is only one of the possible solutions. An equally interesting one is to 
consider a solution heading toward the inflationary attractor a, but still short of it. This would provide closure density 
to a Qm — 0.3 universe, and an acceleration as recently claimed, although at the price of sensibility to the initial 
conditions. Such a model will be investigated in a future work. 
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FIG. 1. Parameter space for the attractors. In the dark-shaded region the solution b is an attractor; in the other regions 
the attractors are the solutions o, c or d, as labelled. In all the paper we focus on the attractor b. The continuous curves mark 
values of Q,/, equal to 0.05,0.1, 0.2, 0.3,...,0.9, top to bottom. The dotted lines are values of p equal to 0.6, 0.65, 2/3, 0.7, 0.8, 
1, 1.3, left to right. There exists also a completely equivalent symmetric region with /3 — » — /3 and /i -^ — /i. 




FIG. 2. Phase space corresponding to fl^ — 0.1, p = 0.7, in MDE. 
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FIG. 3. Contour plot of the exponent m+ of the fluctuation growth law 5c ^ a™ versus n<^,p. The contour levels are for 
m = 0.9 (enclosing the white region), down to in steps of 0.1. In the black region m is complex. Notice that for any given 
il^ the maximum of m is close to p = 2/3, especially for small il^. 
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FIG. 4. Growth of the dark matter fluctuations for various values of the coupling and Q^ = 0.1. In the top panel the trend 
of the horizon length and of two comoving scale (the horizontal lines) show the horizon crossing and the radiation and matter 
eras. Thick lines: wavelength ~ 900Mpc//i. Thin lines: wavelength ~ 100Mpc//i. 
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FIG. 5. Ce spectrum for various models (actually we plot [£{£ + 1)Ci/2it]^''^ , as customary). Notice the shift of the peak 
location for the different values of p, in agreement with the approximation given in the text. The amplitude decreases for 
p 7^ 2/3 (except for values slightly larger than 2/3) and, for a given p, decreases for smaller fl^ , as expected. The data points 
are from Tegmark's home page (http:// www.sns.ias.odu/~max). 
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FIG. 6. Dark matter power spectra A'^(fc) for the same models as above. Notice again that the power is suppressed as long 
as p deviates from the uncoupled law 2/3 (except for values slightly larger than 2/3). Also notice that models with higher p 
are more flattened at small scales. 
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FIG. 7. Contour plot of crsb, ^0), for erg = 1, .9, .8, .7, .6, .5, .4. 
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FIG. 8. Contour plot of x^(p, ^0), for ^ jN = 1, 2, 3, 4, 5, 6. 
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FIG. 9. The parameters in the shaded region have the correct amplitude and shape of the power spectrum, i.e. satisfy the 
constraints ag G (0.45,0.75) and x^/N < 2. 
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